Two-photon speckle as a probe of multi-dimensional entanglement 
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We calculate the statistical distribution P2(/2) of the speckle pattern produced by a photon pair 
current L2 transmitted through a random medium, and compare with the single-photon speckle 
distribution Pi(Ji). We show that the purity Tr of a two-photon density matrix p can be directly 
extracted from the first two moments of Pi and P2. A one-to-one relationship is derived between 
Pi and P2 if the photon pair is in an M-dimensional entangled pure state. For M ^ 1 the single- 
photon speckle disappears, while the two-photon speckle acquires an exponential distribution. The 
exponential distribution transforms into a Gaussian if the quantum entanglement is degraded to a 
classical correlation of M 2> 1 two-photon states. Two-photon speckle can therefore discriminate 
between multi-dimensional quantum and classical correlations. 

PACS numbers: 42.30.Ms, 42.25.Dd, 42.50.Dv, 42.65.Lm 



Optical speckle is the random interference pattern that 
is observed when coherent radiation is passed through a 
difFusor or reflected from a rough surface. It has been 
much studied since the discovery of the laser, because 
the speckle pattern carries information both on the co- 
herence properties of the radiation and on microscopic 
details of the scattering object [1~3]. The superposition 
of partial waves with randomly varying phase and am- 
plitude produces a wide distribution P{I) of intensities / 
around the average (/). For full coherence and complete 
randomization the distribution has the exponential form 
P{I) (X exp(— //(/)). The speckle contrast or visibility. 



V^(/2)/(/)2-l, 



(1) 



equals to unity for the exponential distribution. 

These textbook results [4] refer to single-photon prop- 
erties of the radiation, expressed by an observable /i that 
is quadratic in the field amplitudes. Biphoton optics [5] 
is concerned with observables I2 that are of fourth order 
in the field amplitudes, containing information on the en- 
tanglement of pairs of photons produced by a nonlinear 
optical medium. A variety of biphoton interferometers 
have been studied [6-9], but the statistical properties of 
the biphoton interference pattern produced by a random 
medium remain unknown. It is the purpose of this work 
to provide a theory for such "two-photon speckle" . 

There is a need for a such a theory, because of recent 
developments in the capabilities to produce entangled 
two-photon states of high dimensionality. The familiar [4] 
polarization entangled two-photon state has dimension- 
ality two and encodes a qubit [10]. Multi-dimensionally 
entangled two-photon states include spatial degrees of 
freedom [11-16] and encode a "qudit". The dimension- 
ality of the entanglement is quantified by the Schmidt 
rank M , which counts the number of pairwise correlated, 
orthogonal modes that have appreciable weight in the 
two-photon wave function [17] and is an experimentally 
adjustable parameter [18]. 

As we will show in this paper, two-photon speckle not 
only provides information on the value of M, but it can 
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FIG. 1: Schematic layout (not to scale) of a setup to detect 
two-photon speckle. 



also discriminate between quantum mechanical and clas- 
sical correlations of M modes. For classical correlations, 
on the one hand, the distributions Pi{Ii) and ^2(^2) of 
single-photon and two-photon speckle both tend to nar- 
row Gaussians upon increasing M (with visibilities that 
vanish as \/M). For quantum correlations, on the other 
hand. Pi tends to the same narrow Gaussian while P2 
becomes an exponential distribution. 

We consider a monochromatic two-photon state of elec- 
tromagnetic radiation (density operator pin), scattered 
by a random medium (scattering matrix S). A pair of 
photodetectors in a coincidence circuit is located in the 
far field behind the random medium (see Fig. 1). The 
coincidence detection projects the scattered two-photon 
state (density operator pout) onto quantum numbers k 
and k' , which label the transverse wave vectors of an or- 
thonormal basis of M modes. For a random medium of 
cross-sectional area A and for radiation of wave length A, 
one has N ~ ttA/X^ per polarization degree of freedom. 
The spatial structure of the modes (and the precise value 
of TV) depends on the experimental geometry [17, 18], 
but all our statistical results are independent of it (for 
A/" 2> 1) so we need not specify the modes further for our 
purpose. 

In the far field (at a distance V ^ ^fA from the ran- 
dom medium), the transmitted photon current Ii{k) at a 
given k is detected as a bright spot of area 5 A ~ V"^ /N » 
A^. The random arrangement of bright and dark spots 
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(the speckle pattern) depends sensitively on the realiza- 
tion of the randomness (for example, on the precise con- 
figuration of the scattering centra), and by varying the 
random medium [19] one samples a statistical distribu- 
tion Pi(/i). 

The quantities Ii{k) and Pi{Ii) refer to single-photon 
speckle. The biphoton current l2{k,k') counts the num- 
ber of coincidence detection events per unit time, with 
one photon at k and the other at k' . (We assume k ^ k' .) 
The distribution of I2 in an ensemble of random realiza- 
tions of the disorder is denoted by P2(-^2) and describes 
two-photon speckle. Our goal is to find out what new 
information on the quantum state of the radiation can 
be extracted from P2, over and above what is available 
from Pi. 

The most general two-photon density operator at the 
input has the form 

P-^^ = \HY. ^9i92,gig^4i42|0)(0|Si«9^' (2) 

with oj the photon creation operator in state q and |0) 

the vacuum state. The coefficients in this expansion are 
collected in the N"^ x M"^ Hermitian density matrix p. 
Normalization requires that Trp = 1. If the two-photon 
state is a pure state, then also Trp^ = 1, while more 
generally the purity 

■p = Trp2 g [0,1] (3) 

quantifies how close the state is to a pure state [10] . 

We will present an exact and general theory of the 
speckle statistics for arbitrary pin, and also consider two 
specific simple examples: A maximally entangled pure 
state of Schmidt rank M, 

M 

Ppure = |*m)(*m|, |*m) = M'V^ ^ a\y_^jQ), 

m=l 

(4) 

and its fully mixed counterpart 

M 

Pmixed = «L«-«„|0>(0|«9^«-9..- (5) 

m=l 

Both states (4) and (5) describe a pair of photons with 
anticorrelated transverse wave vectors [20]: If one pho- 
ton has wave vector q^, then the other photon has wave 
vector —qm- (We assume qm ^ for each m.) The dis- 
tinction between the two states is that the two photons 
in state (4) are quantum mechanically entangled, while 
the correlation in state (5) is entirely classical. We will 
see how this difference shows up in the statistics of two- 
photon speckle. 

Scattering by the random medium (in the absence of 
absorption) performs a unitary transformation on the 
creation and annihilation operators. If we collect the 
operators for the incident radiation in the vector a and 
the operators for the scattered radiation in the vector b, 



then b=S-a<^a = S^-b. Substitution into Eq. (2) 
gives the density operator of the outgoing state, 

Pout = i E E P^^<!2,-iW.iS^ ■ • 

91,92 q[,q'^ 

x\0){0\{S^-b),>iS^-b),,. (6) 

From Pout we obtain the biphoton current l2{k, k') by a 
projection, 

hik, k') = la2Tipoutblbl,bkbk', (7) 

where the coefficient a2 accounts for a nonideal detec- 
tion efficiency and also contains the repetition rate of 
the photon pair production. 

We now substitute Eq. (6) into Eq. (7) to arrive at the 
required relation between the biphoton current and the 
scattering matrix, 

l2{k,k') = "2 E E Pliq2,qW-z^kqrSk'q^Sl^,^Sl,^,^. (8) 
91,92 9^,92 

Here we have assumed that p is symmetric in both the 
first and second set of indices, 

/'9l92,9l92 ~ ^92gi,9i92 ~ '^9192,929; ' (^) 

[We can assume this without loss of generality, since any 
antisymmetric contribution to p would drop out of Eq. 

(2).] 

In order to compare with the single-photon current 
Ii{k), we give the corresponding expressions, 

hik) = iaiTVpout&I^fc = "1 E SkqSl,,p^^}, (10) 

9,9' 

in terms of the reduced single-photon density matrix 

^99' = E^992,9'92- (11) 

92 

The coefficient ai is the single-photon detection efficiency 
(which may or may not be different from a2)- 

The next step is to calculate the statistical distri- 
butions Pi, Pj of Ii, l2- Following the framework of 
random-matrix theory [21, 22], we make use of the fact 
that the matrix elements 5*^^ have independent Gaussian 
distributions for 1. (Corrections are of order 1/N.) 
The first moment vanishes, {Skq) = 0, while the second 
moment depends on whether the radiation is detected in 
transmission or in reflection. In transmission through a 
random medium of length L with mean free path I, one 
has 

(15,,]^) = -^^.^. (12) 

Let us begin by calculating the flrst two moments of 
Pi, P2. Carrying out the Gaussian averages, we find for 
the mean values: 

{h) = aia^, {I2) = a2a\ (13) 
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(We omit the arguments k and k' for notational simplic- 
ity.) Neither mean value contains any information on the 
nature of the two-photon state. This is different for the 
variances Var/^ = (if) — (li)'^, for which we find 

Var/i =a?a4Tr(p(i))', (14) 

Var/2 = a^(T^[Trp2 + 2Tr(p(i))^]. (15) 

We conclude that the purity (3) of the two-photon state 
can be obtained from the visibilities Vi = (Var Ji)/(/i)^ 
of the single-photon and two-photon speckle patterns, 



(16) 



This is the first key result of our work. 

To make contact with some of the literature on bipho- 
ton interferometry, we note that in the case of a pure 
two-photon state (when T-" = 1) knowledge of the single- 
photon visibility Vi fixes the two-photon visibility V2. 
The same holds (with some restrictions on the class of 
pure states and with a different definition of visibility) 
for the complementarity relations of Refs. [6-9]. No such 
one-to-one relationship between Vi and V2 exists, how- 
ever, for a mixed two-photon state. 

We next turn to the full probability distribution P2 
of the two-photon speckle. Notice first that, if p is far 
from a pure state, the ratio of the width of the 

distribution and the mean value is ^ 1. Indeed, for 
the fully mixed state (5) one has Trp^j^j,^ — 1/M and 

Tr (pl^^Lod)^ = V2M, so V2 = 2/M « 1 for M > 1. The 
relative magnitude of higher order cumulants is smaller 
by additional factors of 1/M, hence P2 tends to a narrow 
Gaussian for a fully mixed state with M ^ 1. 

The situation is entirely different in the opposite limit 
of a pure state. The density matrix of a pure state fac- 
torizes. 



(17) 



with c a symmetric N'xN' matrix normalized by Tr cc^ — 
1. The corresponding reduced single-photon density ma- 



trix IS p 



(1) 



ccT. The probability distributions P2 and 



Pi in this case of a pure two-photon state are related by 
an integral equation, which we derive in App. A: 



P2(/2) = e(/2) 



Q!2Cr2 



dh ^^^^^^ exp 



ai h 
0.20"^ Ii 



(18) 

Here 6(/) is the unit step function [8(/) = 1 if / > 0, 
6(7) = if / < 0]. 

Without further calculation, we can conclude that 
when Pi is narrowly peaked around the mean (/i), the 
corresponding two-photon speckle distribution is the ex- 
ponential distribution. 



P2ih) oc exp 



ai 



a2(T^ (h 



if Vi < 1. 



(19) 



The limiting exponential form is reached, for example, in 
the pure state (4) for M > 1 (when Vi = 1/2M < 1). 
This is the second key result of our work. 
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FIG. 2: Probability distribution (20) of the two-photon 
speckle for the maximally entangled pure state (4) of Schmidt 
rank M. The exponential distribution (19) (black solid curve) 
is reached in the limit M 00. The black dashed curve shows 
the large-M Gaussian distribution for the fully mixed state 
(5) (plotted for M = 50). 



We can actually give a closed form expression for P2 
in terms of the eigenvalues of the matrix product cc^ 
(see App. B), but it is rather lengthy. A more compact 
expression results for the special case of a maximally en- 
tangled pure state of Schmidt rank M [Eq. (4)] . Then all 
eigenvalues of cc^ are zero except a single 2M-fold degen- 
erate eigenvalue [23] equal to 1/2M. The single-photon 
speckle distribution Pi oc j^^'^-i exp(— 2MJi/q;i(t^) is a 
chi-square distribution with AM degrees of freedom (since 

h T,n=ii\Sk.qJ^ + \Sk.-qJ^) is the sum of 2M Gaus- 
sian complex numbers squared). Substitution into Eq. 
(18) leads to the following distribution of the two-photon 
speckle: 



P2(/2) = e(/2 



X K- 



AM 



2Mh 



a2cr4(2M- 1)! V a2cr4 



Af-l/2 



2M-1 



2Af/2 
a2(T4 



(20) 



The function K2M-1 is a Bessel function. This distribu- 
tion has appeared before in the context of wave propa- 
gation through random media [2] (where it is known as 
the "i^T-distribution" ) , but there the parameter M has a 
classical origin (set by the number of scattering centra) 
— rather than the quantum mechanical origin which it 
has in the present context (being the Schmidt rank of the 
entangled two-photon state). 

We have plotted the distribution (20) for different val- 
ues of M in Fig. 2. The limiting value for /2 ^ equals 



lim ^2(^2) 



2M 



(2M - 1)(72) 



(21) 



The exponential form (19) is reached quickly with in- 
creasing M (black solid curve in Fig. 2). For compari- 
son, we show in the same figure (black dashed curve) the 
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Gaussian distribution reached for large M in the case of 
the fully mixed two-photon state (5). The striking dif- 
ference with the entangled case is the third key result of 
our work. 

In conclusion, we have presented a statistical descrip- 
tion of the biphoton analogue of optical speckle. For an 
arbitrary pure state of two photons, the distribution P2 
of the two-photon speckle is related to the single-photon 
speckle distribution Pi by an integral equation. A nar- 
row Gaussian distribution Pi maps onto a broad expo- 
nential distribution P2. If the two-photon state is not 
pure, there is no one-to-one relationship between Pi and 
P2. For that case we show that knowledge of the visi- 
bilities of the single-photon and two-photon speckle pat- 
terns allows one to measure the purity of the two-photon 
state, thereby discriminating between classical and quan- 
tum correlations of M degrees of freedom. All together, 
this theory provides a framework for the interpretation 
of ongoing experiments [24] on the propagation of multi- 
dimensionally entangled radiation through random me- 
dia. 

We acknowledge discussions with W. H. Peeters and 
J. P. Woerdman. This research was supported by the 
Dutch Science Foundation NWO/FOM. 



APPENDIX A: DERIVATION OF THE 
INTEGRAL RELATION BETWEEN P2 AND Pi 
FOR A PURE STATE 

We start from the expression (8) for the biphoton cur- 
rent, which we write in the more compact form 

l2{k,k') =a2v' ■ M-v'\ (Al) 

in terms of a vector v' with elements v'^ — Sk'q, and a 



we are considering the density matrix p of a pure state, 
which means that M factorizcs as 



matrix M. with elements 



(A2) 
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(We have also defined a vector v with elements Vq = Skq-) 
The probability distribution P2 of I2 is defined by 

P2(72) = {6{l2 -a2v'-M- , (A3) 

with Fourier transform 

/oo 
dl2e'^''P2{l2) 
-00 

= {exp{i^a2v' -M-v'*)),,,'. (A4) 

The average (• • • )yy consists of an average over the scat- 
tering matrix elements Skq contained in the vector v and 
an average over the scattering matrix elements Sk'q con- 
tained in the vector v' . These two averages can be per- 
formed independently, since we have assumed k ^ k'. 

The Gaussian average over v' can be carried out di- 
rectly. 



1 



(A5) 



,Det(l-i^a2CT2A^)// 
To carry out the remaining average over v we need to 
first evaluate the determinant. At this point we use that 



Mqq' = {C-V)qic-V)*,. 



(A6) 



All eigenvalues of this matrix of rank 1 vanish, except 
one nonzero eigenvalue m given by 

f^i=^\{c-v)q\'=v • V* = -Ii{k), (A7) 
9 ^ 

in view of the definition (10) of the single-photon current 
h- 

Because of relation (A7) the average over v in Eq. (A5) 
may be replaced by an average over Ii , 

/•oo 

F2{0 = / dh P{Ii) (1 - i^{a2/ai)a^Ii)~^ . (A8) 
Jo 

Inverse Fourier transformation gives 

n^2) = ^y ^rf^e-«^^P2(0 

= 6(72) 2 / dh ' exp ^— . 

a2(T^ Jo h V "20-2 h ) 

(A9) 

This is the required relation (18) between the distribu- 
tions P2 and Pi of two-photon and single-photon speckle. 



APPENDIX B: CALCULATION OF Pa FOR A 
PURE STATE 

A general expression can be obtained for the distribu- 
tion Pi of single-photon speckle, by repeating the first 
few steps of App. A: 

/oo 
diiS''Pi{h) 
-00 

= (exp(i^Q;i V ■ p^^^ ■ v*))v 
1 



Det(l-i^Q!i<TV(i))' 



(Bl) 



The reduced single-photon density matrix p^^^ has dis- 
tinct positive eigenvalues 7^, each with multiplicity 
and satisfying the sum rule J2ml^mlm = 1- In terms of 
these eigenvalues, Eq. (Bl) can be written as 

m 



(Mm - 1)! 

dXrn 



liXm - i^aia'^jru) ^ 



(B2) 



Inverse Fourier transformation gives the probability dis- 
tribution, 



5 



_1 '\/^Tn-l 



(-1) 



QilO-^Tm (Mm - 1)! 



dXrn 



■ exp 



]J (1 - a;r»7TO'/7m) 



(B3) 



This result for Pi holds for any state of the radiation, but if we now assume that the state is pure, then we can use 
the relation (18) to obtain the two-photon speckle distribution P2 from Pi. (Note that, for a pure state, the "frn's are 
just the eigenvalues of the matrix product cc'.) Substitution of Eq. (B3) into Eq. (18) gives 



a2(T^Jm (Mm - 1)! 



dXrri 



n (1 - a;m7m'/7m)" 



(B4) 



The result (20) given in the main text follows from Eq. (B4) upon taking a single positive eigenvalue 71 = 1/2M 
with multiplicity ii\ = 2M, so that 



P2(72) = 6(72; 



4M (-1) 



2M-1 



a2Cr4 (2M- 1)! 



d 



2M-1 



Kn 2 



2Mxl2 



6(72) 



4M /2Ml2\'^' ^'"^ 



a2a4(2M- 1)! V 0-20"^ 



K. 



2M-1 



2MI2 
0:20-4 



(B5) 
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